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Understanding the cause of the synchronization of population evolution is an important 
issue for ecological improvement. Here we present a Lotka-Volterra-type model driven by 
two correlated environmental noises and show, via theoretical analysis and direct simula- 
tion, that noise correlation can induce a synchronization of the mutualists. The time series 
of mutual species exhibit a chaotic-like fluctuation, which is independent to the noise cor- 
relation, however, the chaotic fluctuation of mutual species ratio decreases with the noise 
correlation. A quantitative parameter defined for characterizing chaotic fluctuation provides 
a good approach to measure when the complete synchronization happens. 
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Synchronization is a common phenomenon in nature (Pikovsky et al. 2001), e.g., in physical 
systems (Pecora & Carroll 1990), chemical systems (Neiman et al. 1999) and biosystems (Winfree 
1990). Specially, the population synchrony has become an important issue in population biology 
(Benton et al. 2001). Additionally, throughout the twentieth century, a dominant study has focused 
on the stochastic fluctuation in population abundance (Cushing et al. 1998; Begon et al. 1996). 
Interestingly, synchronization and stochasticity regarded as two independent phenomena attract 
more and more concerns (Higgins et al. 1997; Leirs et al. 1997; Grenfell et al. 1998; Keeling 2000; 
Tuljapurkar & Haridas 2006). Recently, Benton et al. (2001) suggested that environmental noise 
plays a determined role in the population synchrony. 

Regrettably, due to the difficulties in mathematical analyses and experimental demonstrations, 
the study of relationship between noise correlation and population synchrony is still insufficient. 
However, surprising phenomena emerge sometimes just because there exists correlation or collab- 
oration; one of the most well known examples is the self-organized behavior (Bak et al. 1988). 
We know that two stochastic processes in ecosystems maybe originate from the same source like 
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nature disasters and epidemics. For example, epidemics can lead one kind of species in a mutu- 
alism system as well as another to death synchronously, thus it is reasonable to suppose that the 
population fluctuations of the two species of mutualism system are correlated. 

Mutualism is an association between two species that benefits both. For example, Van der 
Heijden et al. (1998) have argued that mycorrhizal fungi determine plant species diversity. As the 
number of fungal species per system was increased, so the collective biomass of shoots and roots 
increased, as did the diversity of plant community. Namely the number of two species has a positive 
linear relationship, indicating a synchronous fluctuation. For this reason, we focus on the Lotka- 
Volterra-type two-species model with correlated stochastic components. The noise correlation we 
will consider here induces some novel synchronous phenomena in a mutualism ecosystem that are 
not found before. Our main methods will be expected to have applicability to other ecosystems, 
since they have been also of similar correlations of two populations. 

Additionally, noise-induced chaos is observed in many experimental and theoretical situations 
(Miller and Greeberg 1992; Kovanis et al. 1995; Ruiz 1995; Gao et al. 1999; Jing k Yang 2006). 
Can noise induce chaos in a mutualism ecosystem? How does this happen? When is the chaos 
synchronization? In this paper, we will give physical explanations for these problems. 

Mutualistic relationships can be modeled with equations similar to the Lotka-Volterra compe- 
tition equations. 

^ = nx - kix 2 + JiXy, (1) 
dy o 

— = r 2 y- k 2 y + J 2 xy, (2) 

where r y 1 and 7 2 are the positive effects of species y on species x and of species x on species y, 
respectively. r\ and r 2 are respectively growth rate of species x and that of species y. k\, k 2 are the 
parameters related to carrying capacities while one of the mutualists is alone. Due to environmental 
noise like epidemics, natural disasters and adventitious species, realistic mutualism systems do not 
follow the solutions of the above deterministic differential equations exactly. The simplest way to 
include these environmental noises is to add stochastic terms to the above deterministic differential 
equations and the equivalent stochastic differential equations are 

^ = n x- hx 2 + j 1 xy + x£(t), (3) 
-1 = r 2 y- k 2 y 2 +-f 2 xy + yr](t), (4) 

where £(t) and 77 (t) are referred to as the birth rate fluctuations of species x and those of species 
y, respectively. Consider the synchronization of the two fluctuations, we define them as correlated 
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noises. They are Gaussian white noises satisfying 
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in which M x ,M y are the intensities of noises; A, ranges from zero to one, denotes the correlation 
coefficient between and r)(t), i.e., the correlation degree between species x and species y. When 
A equal zero, the growth rate fluctuations of species x and those of species y are independent. When 
A equal 1, their growth rate fluctuations are completely correlated. 5{t — t') is a Dirac delta function 
under different moments. 

In the absence of noises, it can be shown that the system has a stable equilibrium state (xo, 
yo ) when t — > oo for 7 X 7 2 < k\k 2 , where x = (k 2 r\ + r y 1 r 2 )/(kik 2 - J1J2) and Vo = (k\r 2 + 
l 2 f\)lik\k 2 — 7i72)j describes a limit point of system (Stiling 2002). For the sake of simplicity, we 
set k\=k 2 =r\=r 2 =l and r yi='J 2 =0.5. Figure 1 shows that the system can lead to a stable point 
where the isoclines cross. In the presence of noises, we define quantitative parameters, i.e., the 
mean population densities of the mutualists, which are useful to display the influences of noises 
on the mutualism ecosystem. The mean population densities of the mutualists are respectively 
written as 



and their variances are a 2 = (x 2 ) — (x) 2 and a 2 = (y 2 ) — (y) 2 , respectively. Here p(x, y) is the joint 
stationary probability distribution of species x and y satisfying J °° J °° p(x,y)dxdy = 1. 

The equivalent stationary Fokker-Planck equations of Eqs.(3) and (4) can be derived and the 
joint stationary probability distributions can be obtained for two independent noises (Cai and 
Lin 2004). However, they are far difficult to be solved for correlated noises, here we obtain the 
joint stationary probability distributions of the population densities of the mutualists via a direct 
simulation of Eqs.(3) and (4) (Zhong et al. 2006). 

To describe noise-induced chaos and the relationship between chaos and synchrony, we need to 
define chaos carefully. For example, given the time series {x(t)}, an m-dimensional phase portrait 
is reconstructed with delay coordinates, i.e., a point on the attractor is given by {x(t),x(t + 
r), ...,x(t + [m — 1]t)}, with m being the embedding dimension and r being the delay time (here 
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yp(x, y)dxdy, 
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FIG. 1: Graphic model of the mutualism ecosystem is based on modified Lotka-Volterra equations for 
ki—k 2 =r 1 =r 2 =l and 7 1 =7 2 =0.5. There is a stable equilibrium point. 

we choose m = 2, and r = 2) (Wolf et al. 1985; Gao et al. 1999). We denote the initial distance 
between these two points L{to). At a later time t\, the initial length will have evolved to length 
L' (t\). Using Wolf method, we can give lyapunov exponent as following 

A x = (^— Vln-^H) (10) 



in which A is total number of replacement steps. 

It is comprehensible that noises can induce instability. Figure 2(a) illustrates that the popu- 
lations of species x and those of species y driven by environmental noise are instable. However, 
when the noise correlation increases, shown in Figs. 2(b) and 2(c), more and more populations of 
species x and those of species y concentrate on a line in which the ratio of species y to species x 
is constant. When the correlation equals 1, the ratio of species y to species x is independent on 
the noise, i.e., the population fluctuations of species y and those of species x are synchronization. 
This extreme situation can be analyzed from the Eqs.(3) and (4). If the noise correlation equals 1, 
£(i) of Eq.(3)can be regarded as same as i](t) of Eq.(4) for the same noise intensity. If considering 
"# = \%, ~ and substituting Eqs.(3) and (4) for it, we obtain 

dy ^ 

-f = (fci+7 2 )j/-(*2 + 7i)- (11) 
at x 

Considered a steady solution -J- = 0, then Eq(ll) can be written as | = = 1, i.e., the 

populations of species x have a positive linear relationship with those of species y, shown in 

Fig.2(c). 
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FIG. 2: Distributions of the mutualists driven by environmental noise with different correlations, (a) A = 0.0; 
(b) A = 0.8; (c) A = 1.0, including simulated results (circle) and theoretical analysis (solid line) results. The 
parameter are M x = M y = 0.8. 

The linear relationships between the variance of | and the noise correlation are given in Fig. 3, 
indicating the variance of | decreases with the noise correlation. When the correlation equals 1, 
the fluctuations of the ratio of species y to species x are extinct. 

Figure 4 shows that the lyapunov exponents of species y and those of species x fluctuate between 
0.2 and 0.35 when the noise correlation increases, however, the lyapunov exponent of the ratio of 
species y to species x drops to zero. Generally speaking, when lyapunov exponent is positive, the 
time series are regarded as disorder or chaotic characteristic. On the contrary, if lyapunov exponent 
is non-positive, the time series are of order or periodic characteristic. Figure 4 indicates that two 
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FIG. 3: Variances of the ratio of species y and species x vs. noise correlations under different noise intensities. 
The remaining parameter are the same as for Fig. 2. 



FIG. 4: Lyapunov exponents of species y, species x and their ratio vary with noise correlations. The 
completely synchronous fluctuation point corresponds to zero-lyapunov exponent. 

disorder time series with correlation can exhibit an order behavior, and the synchronization in a 
mutualism ecosystem happens at the lyapunov exponent equals zero. 

In summary, the environmental noise correlation can cause a synchronous fluctuation in mutual- 
ism ecosystems, and this synchronous fluctuation can be characterized by the lyapunov exponent. 
If environmental noise correlation originates from the correlation degree of the mutualists, our 
main results illustrate that the mutualism ecosystems with closer relationship have more stable 
evolutions. This is a helpful illumination for ecologists to study ecosystem balance. 
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